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1. Introduction 

Star formation is one of the key processes in the Universe and is of fundamental importance to 
astrophysics (McKee & Ostriker 2007). We still do not have a full understanding of star formation 
and, in particular, on the role that magnetic fields play (Crutcher 2012; Li et al. 2014). Cosmic Rays 
(CRs) are the main ionising and heating agent in dense, starless, molecular cloud cores (Padovani 
& Galli 2013). The role of magnetic fields has been debated for decades, but due to a lack of 
precise measurements it has still not been conclusively settled. There arc two major classes of 
star-formation theories: i) strong field models where the B-field controls the molecular cloud, with 
ambipolar diffusion driving the formation and collapse of dense cores, and ii) weak-field models 
where turbulent flows trigger star formation. Although B-fields arc at least comparable in strength 
to the turbulent pressure in the ISM, there is no definitive evidence for B-fields dominating gravity 
or for ambipolar-diffusion-driven star formation (Crutcher 2012). 

It is very difficult to measure the detailed properties of dense molecular clouds, and perhaps 
the most difficult property is the magnetic field - its strength and geometry. The most common 
method is via Zeeman splitting of radio-frequency HI, OH, and CN lines or masers, which typi¬ 
cally gives only the magnetic field strength along the line-of-sight, /i| os (Crutcher 2012) 1 . Other 
methods arc also difficult, which include using measurements of optical and near-infrared polar¬ 
isation (extinction along the line-of-sight) and sub-mm polarised thermal dust emission (difficult 
from the ground) (Poidevin et al. 2013). 

An alternative method of probing the magnetic field is via synchrotron radiation (Brown & 
Marschcr 1977; Marscher & Brown 1978; Orlando & Strong 2013), which has rarely been men¬ 
tioned in the literature in relation to molecular clouds. Synchrotron radiation is produced primarily 
by relativistic CR electrons when decelerated by magnetic fields. The intensity of synchrotron ra¬ 
diation depends only on the number and energy spectrum of CR electrons and the magnetic field 
strength perpendicular to the line-of-sight. We know the measured energy spectrum of CRs (Ack- 
ermann et al. 2012), at least on average at energies of relevance to radio synchrotron emission 
(~GeV), and we know that the CR density varies slowly throughout the Galaxy and can pene¬ 
trate dense molecular clouds at these energies and above (Brown & Marschcr 1977; Marschcr & 
Brown 1978; Umebayashi & Nakano 1981). 2 In molecular clouds, the magnetic field strength is 
much larger than in the ambient interstellar medium (Crutcher 1999). Therefore, in principle, the 
synchrotron intensity should give a detectable signature, which could be used as a probe of the 
magnetic field. This is a direct way of measuring the total magnetic field strength, including the 
irregular (turbulent) component, which most other indicators (Zeeman, optical polarisation) are 
not directly sensitive to, since they measure the regular (Zeeman) or ordered (optical polarisation) 
field. 3 Polarised synchrotron emission could provide additional information, including the ordered 

1 It is possible to get the total B-field strength when complete line splitting is observed, which is possible with masers 
(Crutcher 1999); see also Robishaw (2014). 

2 There have been claims that GeV CRs cannot fully penetrate the densest clouds, thus suppressing the CR diffusion 
coefficient (Jones et al. 2008; Protheroe et al. 2008; Jones et al. 2011). 

3 Other indirect measures of the random component exist, such as the Chandrasekhar-Fermi method, which uses 
the dispersion of the measured polarisation angles to probe the magnetic field in the plane of the sky (Chandrasekhar 
& Fermi 1953; Watson et al. 2001; Crutcher et al. 2004). See also Hildebrand et al. (2009) and references therein for 
further extensions. 
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versus anisotropic random component and projected angle of the magnetic field on the sky. 

It is therefore somewhat surprising that very little attention has been given to using synchrotron 
as a probe of molecular clouds. Jones et al. (2008) observed two nearby (3-4 kpc) dense cold 
starless cores (G333.125-0.562 and IRAS 15596-5301) with the Australian Telescope Compact 
Array (ATCA) at 1384 and 2368 MHz, to try to detect secondary leptons. 4 They found upper 
limits of ~ 0.5 mJy/beam and constrained the B-field strength to B < 500 /./G. However, this is still 
compatible with the scaling of |B| and «h - more sensitivity is required. Protheroe et al. (2008); 
Jones et al. (2011) only found upper limits from SgrB2 after subtraction of the dominant thermal 
emission. The only possible candidate so far is from the GO. 13-0.13 molecular cloud detection, 
which was detected at 74 MHz, with an associated CO hotspot (Yusef-Zadeh et al. 2013). However, 
a displacement between the radio position and molecular core suggests it could be from a different 
region of space. 

These non-detections can be partly understood due to the relatively weak (typically mJy or 
less) signal that is expected to come from in-situ synchrotron emission inside the cloud itself. 
This is due to the fact that on large scales (~ 1-10 pc), the magnetic fields in clouds appeal - to 
be relatively weak (~ 10 /./G) while strong fields are on scales much smaller than this (~ 0.05 pc) 
resulting in a weak flux signal. Also, most of the collapsing clouds ("cores") are located at low 
latitudes where there is significant confusion from background synchrotron and free-free emission. 
Nevertheless, high resolution and high sensitivity observations could allow molecular clouds to be 
mapped in some sight-lines. This may also shed light on CR penetration into the densest clouds, 
which sometimes appear as a radio dark cloud (RDC) (Yusef-Zadeh 2012). Note that recent high- 
resolution 5/20 GHz JVLA observations of the Galactic centre cloud G0.216+0.016 have detected 
compact (< 2.2 arcsec, or sub-pc) non-thermal sources, which may be the signature of in-situ syn¬ 
chrotron radiation from secondary CR electrons (Jones 2014). 

In this chapter we briefly review the physics of synchrotron radiation and magnetic fields, and 
the relation that appeal's to exist between them in molecular clouds. We then give the prospects of 
detecting this signature in molecular clouds with the SKA. 

2. Synchrotron radiation 

Synchrotron radiation is emitted primarily by relativistic cosmic-ray electrons spiralling in the 
Galactic magnetic field. It is this radiation that often dominates the radio sky at frequencies below 
a few GHz. The theory of synchrotron radiation is well understood. For a power-law distribution 
of electron energies, 


N(E)dE =NoE~' y dE , (2.1) 

the emissivity, j v , of synchrotron radiation is given by 

j v « A 0 B (7+1)/2 V (1 “ r ) /2 , (2.2) 

4 In this article we focus on primary CR electrons, although the conversion into secondary leptons could be signif¬ 
icant in some clouds (Dogel’ & Sharov 1990; Protheroe et al. 2008). Observations with the SKA could also provide 
useful constraints on secondaries. 
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where B is the magnetic field strength, /Vo is the number density of electrons, y is the power-law 
index of electron energies, and v is the observing frequency. Highly relativistic (GeV and above) 
CR electrons are expected to penetrate dense clouds freely (Umebayashi & Nakano 1981). If 
the electron energy spectrum is a power-law with slope — y then the observed synchrotron radio 
emission spectrum is also a power-law with slope a = (1 — y)/2 (flux density S °c v“). The cosmic 
ray energy distribution at energies of order GeV can be approximated by a power-law with slope 
y« +2.5-3.0 (Ackermann et al. 2012), which corresponds to a synchrotron index a « —0.8. This 
is indeed the typical spectral index observed at GHz frequencies (Reich & Reich 1988; Platania 
et al. 1998). 

It can also be seen that the emissivity scales as r 0'+ 1 )/ 2 , which means it goes as approxi¬ 
mately B 2 . This is of relevance to molecular cloud collapse, since the magnetic field is expected 
to be significantly amplified during collapse, and thus could give a detectable signal from in-situ 
synchrotron radiation. 

3. Magnetic fields in collapsing clouds 

The role that magnetic fields play in the processes of molecular cloud and star formation has 
been debated for decades. Theoretical studies suggest that magnetic fields play an important if 
not crucial role in the evolution of interstellar clouds and the formation of stars. In summary, 
magnetic fields provide magnetic support against cloud collapse. There are various models of star 
formation and the details of the magnetic field arc always important. For example, the core of a 
cloud can become unstable due to ambipolar diffusion, collapsing to form stars, while the envelope 
can remain in place. The connection between the core and the surrounding envelope by magnetic 
field lines can transfer angular momentum outward and make it possible for stars to form. Other 
star formation models have the dissipation of magnetised turbulence as a controlling factor in star 
formation. Measuring the magnetic field is a key observation that allows us to infer i) whether 
supersonic motions are Alfvenic, and ii) the relative importance of the gravitational, kinetic and 
magnetic densities in dense clouds (Crutcher 1999). These observables thus allow us to test star 
formation models such as ambipolar diffusion and turbulence (Crutcher 2012; Lazarian et al. 2012). 

Detailed measurements of the magnetic field strength and alignment are difficult. However, 
in recent years, direct measurements of the magnetic field strength have been made. Most notable 
are Zeeman splitting data (Crutcher 1999; Crutcher et al. 2010) and also sub-mm thermal dust 
emission (Poidevin et al. 2013). Detailed studies of Zeeman splitting from a sample of molecular 
clouds indicate that the thermal-to-magnetic pressure f5 p ~ 0.04, implying that magnetic fields 
are important. Moreover, the measurements showed that magnetic field strengths scale with gas 
densities as B °< n K « 0.5—0.7, as shown in Fig. 1. This is close to the theoretical value K = 
0.47 predicted by models of ambipolar diffusion (Fiedler & Mouschovias 1993). The latest value 
appeal's to be K = 0.65 (Crutcher 2012) but there is considerable scatter in the measurement (Fig. 1); 
our best-fitting value applied to detections above 3a is re = 0.54 + 0.05, although there could be 
biases when neglecting non-detections (Crutcher et al. 2010). The large scatter may be related to 
the fact that Zeeman splitting is only sensitive to the regular (ordered and directional) magnetic 
field component along the line-of-sight; the B-«h relation may be different for turbulent fields. 
Furthermore, this trend only occurs above some density «o ~ 300 cm 3 , although this has still to 
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be determined precisely. Clearly more data, and complementary probes of the magnetic field, are 
needed to make progress in this area. 



Figure 1: Zeeman splitting measurements of the magnetic field of a number of molecular clouds, plotted 
against the volume density of molecular gas. Significant (> 3c) detections are shown as black filled circles 
(data taken from Crutcher et al. 2010). There is considerable scatter, yielding various slopes (see overplotted 
lines) depending on the exact model being fitted. The power-law slope between B z and density nn is in the 
range K « 0.5-0.7. Our best fit (solid line) for detections greater than 3 <J significance yields K = 0.54 ± 0.05 
above «n = 300 cm 3 . Preliminary figure reproduced from a forthcoming publication (Strong et al. 2015). 


4. Predictions for synchrotron radiation from collapsing clouds 

Given that the synchrotron emissivity scales as ~ B 2 and B scales as ~ «^ 6 , it is logical that 
it should also scale roughly as the volume density i.e. j v <=< /7 H . From this, one might expect 
low frequency maps such as the Haslam et al. 408 MHz map (Haslam et al. 1982) to be bright 
around giant molecular clouds (GMCs) and for molecular clouds to be very bright in high resolution 
observations (e.g. VLA, ATCA). We will now use the observed scaling relation of B with /?h to 
estimate the flux density expected for typical molecular clouds. 

We assume that the ambient CR electrons pervade molecular clouds unimpeded and a power- 
law distribution of CR electron energies with slope y and a power-law slope between density nu and 
B-field strength B above a value no = 300cm 3 . Using the CR flux model of Strong et al. (2011), 
the predicted brightness temperature (in mK) at 408 MHz can be be approximated by (Strong et al. 
2015): 


/ Tpred 

mK 


= 2.8 x 10 3 


Nu 

10 23 cnr 


» H 

300 cm 


X K(r+l)/2-l 


(4.1) 


where Ah is the column density (cm 2 ) and «h the volume density (cm 3 ). This corresponds to a 
predicted integrated flux density (in mJy) at 1 GHz, for a source subtending a solid angle n src , 
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Table 1: Molecular cloud data from Crutcher (1999) with estimates of predicted integrated synchrotron flux 
density .S'giiz (mJy) at 1 GHz based on the statistical B-nu scaling law (equation 4.2) with K — 0.6. The flux 
density has been scaled to 1 GHz assuming a spectral index a = —0.8 ( 7 = 2.6). The brightness temperature 
7 qhz is what would be observed with a 51 arcmin beam. 


Name 

[pG] 

nu 2 

[cm -3 ] 

R 

[pc] 

D 

[kpc] 

6 

[arcsec] 

7ghz 

[mK] 

SgHz 

[mJy] 

W3 OH 

3100 

6.31 x 10 6 

0.02 

2.0 

4.0 

0.06 

0.05 

DR210H1 

710 

2.00 x 10 6 

0.05 

1.8 

11.2 

0.30 

0.27 

Sgr B2 

480 

2.51 x 10 3 

22.0 

7.9 

1149 

1200 

1000 

M17SW 

450 

3.16 x 10 4 

1.0 

1.8 

236 

31 

27.0 

W3 (main) 

400 

3.16 x 10 5 

0.12 

2.0 

24.3 

0.49 

0.43 

S106 

400 

2.00 x 10 5 

0.07 

0.6 

48.1 

0.74 

0.65 

DR210H2 

360 

1.00 x 10 6 

0.05 

1.8 

11.2 

0.14 

0.13 

OMC-1 

360 

7.94 x 10 5 

0.05 

0.4 

50.3 

2.3 

2.0 

NGC2024 

87 

1.00 x 10 5 

0.2 

0.4 

196 

14.6 

13.0 

W40 

14 

5.01 x 10 2 

0.05 

0.6 

34.4 

0.04 

0.03 

p Oph 1 

10 

1.58 x 10 4 

0.03 

0.1 

91.7 

0.14 

0.13 



= 6.6 x 10 6 



/ Nu \ / »h \ fc(y+l)/2-l 
\ 10 23 cm -2 / \300cm~ 3 / 


(4.2) 


Table 1 lists some example molecular clouds, using data from Crutcher (1999), with predicted 
flux densities at 1 GHz. We have used the B — nu relation above with K = 0.6, assume a synchrotron 
frequency spectral index a = 1.0 ( 7 = +3.0), and G src = 71 / 4 x 0 2 . Dense molecular clouds have 
typical densities of 10 5 -10 6 cnr 3 in H 2 and linear sizes of ~ 0.05 pc. This gives column densities 
of ~ 10 23 cm -2 . For typical distances of a ~kpc, this corresponds to angular sizes of ~ 10 arcsec. 
It can be seen that many of these sources have predicted flux densities of ~mJy. It is interesting 
to see that a few sources have much larger predicted flux densities (e.g. SgrB2 at about 1 Jy). 
However, one has to be careful since these are due to the large physical size assumed (22 pc for 
SgrB2). In practice, the collapsing clouds tend to be very small, often clustered, in a parent cloud 
that is much larger. The magnetic field measured in the densest regions is unlikely to apply to the 
entire cloud. Thus it is easy to over-estimate the flux density in this way and this appears to be 
why dense molecular clouds are not bright in low resolution radio surveys such as the Haslam et al. 
(1982) 408 MHz map. On the other hand, additional synchrotron from secondary leptons could 
boost the synchrotron level (Protheroe et al. 2008). 

Therefore, the predicted flux densities should only be considered order-of-magnitude estimates 
at this point since the precise values depend very sensitively on the choice of fc and no and on the 
observed input parameters. Furthermore, the huge scatter about this relation observed in Fig. 1 al¬ 
ready indicates that either the measurements are not representative of the mean field, or, the simple 
B — nu relationship does not hold. New observations will be crucial for testing this hypothesis. 
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5. Prospects for the SKA 

The high sensitivity and high resolution of the SKA should finally provide the first definitive 
detection of in-situ synchrotron radiation from molecular clouds themselves. The dense cores 
inside clouds are typically 0.05 pc in size, corresponding to ~ 1 arcsec at a distance of ~kpc. 
This is well matched to the core of SKA1-MID, which provides good u,v coverage (and therefore 
good surface brightness sensitivity) up to ~ 5 arcsec. Estimated flux densities are of order mJy for 
typical dense molecular clouds at GHz frequencies (see Table 1), which is approximately the depth 
to which current observations have been targeted (e.g. Jones et al. 2008). Optimal frequencies 
are in the range of a few hundred MHz (band 1-2 of SKA1-MID) to ~ 2 GHz (band 3 of SKA1- 
MID). The very lowest frequencies (< 300 MHz) will be affected by synchrotron and free-free 
self-absorption. Above a few GHz, the steep spectrum (a « —0.8) means the signal will be very 
weak and more susceptible to free-free contamination, which has a flatter spectral index (a & —0.1 
in the optically thin regime). The SKA is therefore the ideal instrument for studying this, yet 
untapped, area of radio astronomy. 

Multi-frequency (matched u,v coverage) continuum mapping will be important to verify that 
the emission is non-thermal. In principle, optically thin thermal emission could be removed using 
radio recombination lines, which could be detected using the wide-band spectral capability of the 
SKA. Mapping the polarisation of the synchrotron radiation provides further information on the 
magnetic field, such as the orientation and orderliness. The polarisation fraction could be up to 
~ 75 % for a ordered field with no depolarisation. However, in practice the observed level will be 
lower, making this a more difficult measurement to make. Lower frequencies (bands 1/2) will also 
likely be affected by Faraday depolarisation along the line-of-sight. 

Diffuse clouds (e.g., Heiles 1997; Myers et al. 1995) will often not be bound structures, have 
modest magnetic fields and an in-situ synchrotron signal will not be obvious due to confusion from 
nearby objects. Given their large size (arcmin to several degrees) and diffuse nature, they will not 
provide the cleanest signal and are not well-suited to SKA baselines. More well-defined clouds 
such as CO/sub-mm catalogues of molecular clouds (e.g., Roman-Duval et al. 2010; Planck Col¬ 
laboration et al. 2011) provide lists of hundreds of clouds but, due to the relatively low angular 
resolution, they are relatively large in size (5-20 pc) and may not contain the densest clouds of in¬ 
terest here. Proto-stellar outflows may also emit appreciable synchrotron radiation (Sokoloski et al. 
2008) which will need to be carefully mapped. Nevertheless, with the large field-of-view of the 
SKA, it should be possible to find useful areas (Galactic plane and molecular clouds, particularly 
at high Galactic latitude where cleaner sight-lines should exist) to survey to look for the signature 
of in-situ synchrotron emission. 

The best targets will be the most compact and dense clouds, where the magnetic fields are 
likely to be highest, and the signal will be the cleanest. Examples such as those in Table 1 , with 
angular sizes of ~arcsec are ideal for the SKA. An example would be the ring of dense molecular 
cores in the W40 complex, where upper limits of ~mJy exist but require deeper observations 
(Pirogov et al. 2013). The long baselines of SKA will in fact allow us to resolve the emission 
inside a dense cloud opening up the potential of studying the magnetic field structure within the 
cloud. Surveying a number of clouds in this way is likely to directly shed light on models of star 
formation through the alignment of the magnetic field. They should ideally have as little or no 
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free-free emission from ionised gas or nearby supernova remnants (SNRs). 5 In many cases, such 
surveys could piggy-back on other SKA surveys of the Galaxy. 

The high density and quiescent nature of infrared dark clouds (e.g., Peretto & Fuller 2009), 
together with high resolution column density images available, will make them interesting targets 
for SKA observations to probe the magnetic field at early stages of formation. A good example 
would be the massive star-forming cores within the infrared dark cloud SDC335.579-0.272 (Peretto 
et al. 2013). 

6. Conclusions and outlook 

Observations of the properties of dense molecular clouds arc critical in understanding the 
process of star-formation. One of the most important, but least understood, is the role of the 
magnetic fields. We propose to use high resolution, high sensitivity radio observations with the 
SKA to measure the in-situ synchrotron radiation from these molecular clouds, complementing 
other methods such as Zeeman splitting (Robishaw 2014). If the CR particles penetrate as expected, 
then we can measure the B-field strength directly using radio data (the CR flux is relatively smooth 
and varies slowly throughout the Galaxy with a scale-height of ~ 1 kpc; Orlando & Strong 2013; 
Stepanov et al. 2014). If they cannot, then the ability of the SKA to pick up extended, low surface- 
brightness emission may allow us to search for synchrotron emission from the outer regions of such 
clouds. Collapsing cores are typically arcsec in size and flux density estimates are ~ mJy at 1 GHz 
and thus should be readily detectable by the SKA. The large field-of-view will allow many lines-of- 
sight to be investigated for a given molecular cloud. Multiple frequencies and radio recombination 
lines will allow separation from free-free emission while polarised data opens up the possibility of 
mapping the field geometry within collapsing clouds. 

CD acknowledges support from an ERC Starting (Consolidator) Grant (no. 307209). 
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